Evidence of Cooper pair pumping with combined flux and voltage control 
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We have experimentally demonstrated pumping of Cooper pairs in a single-island mesoscopic 
structure. The island was connected to leads through SQUID (Superconducting Quantum Inter- 
ference Device) loops. Synchronized flux and voltage signals were applied whereby the Josephson 
energies of the SQUIDs and the gate charge were tuned adiabatically. From the current-voltage 
characteristics one can see that the pumped current increases in le steps which is due to quasipar- 
ticle poisoning on the measurement time scale, but we argue that the transport of charge is due to 
Cooper pairs. 
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A device that yields a DC current in response to an AC 
signal at frequency / according to the relation / = Qf 
is called a charge pump. In the case of electron pumps 
Q = me while for Cooper pair pumps Q = 2me, where 
m is an integer denoting the number of charges being 
pumped per cycle. Typically pumping electrons in meso- 
scopic structures requires an array of at least three tun- 
nel junctions with voltage gates coupled to the islands 
in between the junctions. A Cooper pair pump is ob- 
tained when the tunnel junctions are replaced by Joseph- 
son junctions. These devices appear at first sight to be 
very similar and actually the very same samples may 
serve as both Cooper pair and electron pumps depend- 
ing on whether the device is in the superconducting state 
or not. However, major differences exist. Besides the 
doubled charge in the superconducting state, the nature 
of the tunneling processes is very different, too. Elec- 
trons can tunnel downhill in energy due to the inherent 
dissipation mechanisms in normal metals with the rele- 
vant time scale given by the RC time constant, where 
R is the tunnel resistance and C the tunnel capacitance. 
Cooper pairs, on the other hand, try to conserve their 
energy, and in the absence of an electromagnetic envi- 
ronment, (i.e. zero impedance) only elastic processes are 
possible. Their maximum pumping frequency is propor- 
tional to Ej/(EcK), where Ej and Eq are the Joseph- 
son and charging energies, respectively. What is more, 
superconducting circuits may behave coherently in the 
quantum-mechanical sense. The first attempt to pump 
Cooper pairs dates back to over a decade agoi However, 
Cooper pair pumps have not been even nearly as accu- 
rate as single-electron pumps. The best example of the 
latter ones is the NIST seven-junction pump 2 . The mo- 
tivation behind pumping Cooper pairs is two-fold. First 
of all, Cooper pair pumps are hoped to be able to pump 
larger currents than their normal state counterparts while 
still being accurate. This is roughly because increasing 
E"j/(EcK) is easier than increasing 1/{RC). Secondly, 
the operation of Cooper pair pumps is interesting from 
the point of view of secondary "macroscopic" quantum 
phenomena and the structures are quite similar to the 
superconducting qubits (see, e.g., Refs-HQ). Pumping of 



electrons using surface acoustic waves is another active 
field of study, see, e.g., Ref. @- 

In this work we report on the experimental demonstra- 
tion of pumping Cooper pairs in a structure nicknamed 
the Cooper pair "sluice" introduced and theoretically an- 
alyzed recently by us, see Ref. 0. The device is particu- 
larly simple; it has just one superconducting island, like 
the single Cooper pair transistor, but the bare Josephson 
junctions are replaced by SQUID loops. The device may 
be alternatively viewed as a tunable Cooper pair box, a 
Josephson charge qubiti Here the control is achieved 
via adiabatically manipulating both the fluxes through 
the two loops and the gate voltage. Ideally the SQUIDs 
act as tunable Josephson junctions whose coupling en- 
ergy can be varied between a value close to zero and the 
sum of the couplings of the individual junctions. First we 
describe the experimental setup and discuss the theoret- 
ical idea briefly. Then we present measured data of the 
pumping experiment. We demonstrate that the pumped 
current obeys nicely the theoretical predictions. We also 
comment on possible ways of improving the results should 
the device be used in applications and discuss the signif- 
icance of the results. 

Figure ^ shows an SEM image of the sample used in 
the experiments along with a schematic of the measure- 
ment setup in Fig^c). The device was fabricated out of 
aluminum using standard e-beam lithography and two- 
angle shadow evaporation. It consists of a superconduct- 
ing island that connects to the leads via SQUID loops. 
These are relatively large (10 fim by 100 fiva) in order 
to have good inductive coupling but the island and the 
junctions are still small such that the charging energy is 
large enough («I K) to suppress thermal effects. The 
sample was attached to a dilution cryostat with a base 
temperature of 20 mK with the RF-lines connected. 

Ideally, the pumping of m Cooper pairs is achieved by 
applying the three pulses in Fig[2^b) through the atten- 
uated RF-lines. The upmost signal is applied to the gate 
while the two lower ones represent the currents flowing in 
the input coils. Two different versions of the gate pulse 
are shown, one for pumping "forward" and one for pump- 
ing "backward" . To understand how the device works, 
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FIG. 1: (a) Scanning electron micrograph of the sample. The 
two input coils can be seen on top and bottom, respectively. 
The gate extends to the far right and the gate capacitance is 
C g = 0.24 fF based on DC measurements. The current flows 
between the two leads on the left side, (b) Closeup of the 
island. The measured total capacitance of the island is 3.7 
fF which corresponds to a charging energy of about 1 K for 
Cooper pairs. The maximum Ej per SQUID is estimated to 
be around 0.5 K based on the normal state resistance, (c) 
Schematic illustration of the measurement setup. We used 
commercial room temperature electronics for the current mea- 
surement and three synchronized arbitrary waveform gener- 
ators for the control pulse. The external coil for tuning the 
background of the SQUIDs is at 20 mK. The voltage biasing 
happens via voltage division through resistive lines. A surface 
mount capacitor of 680 pF and an on-chip capacitor on the 
order of 10 pF were also used. 



it is instructive to look at the Hamiltonian of the device, 
which reads 



H =E c (n - n g ) 2 ~ cos(0 + <p/2) 

-£?($ 2 )cos(p/2-$. 



(1) 



Here Eq — 2e 2 /Cs is the charging energy for Cooper 
pairs where Cs is the total capacitance seen from the is- 
land. Furthermore, E\ with j = 1,2 are the (signed) 
Josephson energies of the two SQUIDs which can be 
tuned with the external fluxes $j. For identical junctions 
E\ = Ef a * cos( 7 r$ J /$ ), where $ w 2 x lfr 15 Wb is the 
flux quantum and E™ ax is proportional to the critical cur- 
rent Jc of the individual junctions via Ef aK = (K/e)Ic- 
Furthermore, n g — C g V g /2e is the gate charge in 2e units, 
h is the number operator for Cooper pairs, <p is the phase 
on the island and their commutator is [h, <f>] — i. The en- 
vironment couples to the pump through ip which is the 
phase difference over the pump. If the SQUIDs were 
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FIG. 2: (a) Contour plot of the measured DC current at con- 
stant voltage against DC currents in the two input coils. The 
total variation in the current is around 40 pA at this bias point 
(150 jiiV). The arrow line indicates the path along which the 
flux pulsing is performed in the pumping experiment. The 
lines of minimum current along which the arrows are aligned 
are the lines along which half a flux quantum threads one of 
the two SQUIDs. The slight tilting of the lines is a signa- 
ture of the inductive cross-coupling. Arranging the pulsing 
as shown compensates for the cross coupling, (b) Waveforms 
that were used in the experiment. The thin almost sinusoidal 
pulse is the gate signal for pumping in, say, "forward" direc- 
tion, and the dashed 7r-shifted signal is for pumping in the 
"backward" direction. The low level of the gate pulse is zero. 
The thick lines are the current signals corresponding to the 
arrowed path in the previous contour plot, (c) Contour plot of 
the measured current at a constant voltage of 250 fiV against 
the relative phase differences between the signals with the 
pumping signal being applied at 2 MHz. The blue circle is 
the optimal choice for pumping "forward" while the red circle 
is the optimal point for pumping "backward". The ampli- 
tude was set large (over 400e) and the variation in current 
was 150 pA. This operation point is far from optimal, but we 
still obtain a clear modulation for calibration purposes. 



to have perfectly identical junctions as well as vanishing 
sclf-inductance and if the flux control were perfect then 
the effective Josephson couplings could be set to zero. 

Figure[2Ia) shows a contour plot based on the measure- 
ment of the current through the device at a constant volt- 
age against the DC currents in the two input coils. Along 
the lines of minimum current the flux through either of 
the loops is (k + l/2)$ , where k is an integer. The mea- 
surement reveals not only the mutual inductances My 
between coil i and SQUID j, which were Mn = 30 pH, 
M12 = 2 pH, M 2 \ = 3 pH and M 22 = 50 pH, but also 
the proper offsets at any given time, i.e. the background 
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fluxes threading the loops. This measurement does not 
fully demonstrate to which extent it is possible to sup- 
press the Josephson energy. 

In the beginning of an ideal pumping cycle the Ej's 
of both loops are set as close to zero as possible and the 
position of the gate determines the ground state. We see 
that initially the ground state of the island is an eigen- 
state of charge. We then adiabatically "open" one of the 
SQUIDs, i.e. move to the tip of the, say, horizontal ar- 
row in Fig E{a) which means that the Ej of the SQUID 
1 is maximized while for the other it is still zero. We 
stay at the tip of the arrow for some time and start to 
either decrease or increase the gate charge n g depending 
on the direction we have chosen. When the gate reaches 
its extremum we "close" the SQUID again. Now if every- 
thing has been adiabatic the system is still in its ground 
state. The charge is again a good quantum number at 
this point but since the position of the gate is different, 
the number of charges is different too. The only possi- 
bility is that the excess charges have tunneled through 
the SQUID whose Ej has been non-vanishing during the 
cycle. The Ej of the second SQUID is then opened and 
the gate put back to its initial position. Finally the sec- 
ond SQUID is also closed. The number of Cooper pairs 
pumped is given by the difference between the integers 
closest to the high and low level of the gate charge. Fix- 
ing the low level and sweeping the high level should result 
in a 2e-periodic staircase in the pumped current. 

The phase of the gate determines naturally the direc- 
tion, i.e. a 180-degree phase shift reverses the pumped 
current. Fig^c) illustrates the measured behavior of the 
current when the relative phases between the pulses are 
varied. The phase of coil 1 is fixed at 180 degrees and 
the phases of the other two are swept. The two circles 
shown are the optimal choices for pumping. Note that 
the extrema of current are indeed 180 degrees apart in 
the gate as expected and the optimal choices are the ones 
illustrated in Fig. EI F° r practical reasons we were forced 
to use frequencies in the MHz range, but in the present 
pumping scheme it is possible to increase the value of 
current conveniently by increasing the gate amplitude. 
We tried out different shapes of pulses such as a mere 
sinusoidal gate signal, but it was found that it is better 
to keep the gate constant while the Ej is not maximized 
which is in accordance with the adiabaticity requirement. 
In practice we have arranged for a 15% dead time between 
the flux pulses although no systematic optimization of 
the pulses was performed. 

Figure Ola) shows an example of characteristic IV- 
curves (i.e., current-voltage curves) with the pumping 
signal being applied at / = 3 MHz. The effect of the 
change of direction is shown. The curves correspond to 
eight different values of gate amplitude. We see immedi- 
ately that a leakage current exists on top of the pumped 
current that is on the same order or less than the pumped 
current. The IV-curves, however, clearly shift and the 
curves for pumping in opposite directions are far apart. 
The total current flowing through the device is a sum 




FIG. 3: (a) Examples of measured IV-curves with the pump- 
ing signal applied at 3 MHz. The gate charge (in 2e units) 
varies between 4 and 34. The solid curves correspond to 
pumping forward and the dashed curves correspond to pump- 
ing backward. Here V mca s is the measured value of voltage 
over the pump, (b) Difference of current, AI, in the IV- 
curves of (a) for pumping in opposite directions. The dotted 
lines indicate the expected values. 



of two contributions, one being the leakage supercurrent 
that can be associated with the dynamical phase of the 
wave function and the other being the less trivial pump- 
ing contribution attributable to the geometric phase. If 
one assumes that the leakage is the same for the pumping 
in both directions at a definite voltage bias point, then 
the difference between the IV-curves should be twice the 
magnitude of current pumped in this case. Fig. re- 
veals that at low voltages (tens of //V) and at smaller 
amplitudes this pumping contribution is indeed close to 
the expected level shown with dotted lines. The leak- 
age current which is due to the nonideal environment 
and flux control is undesirable from an application point 
of view, but the physical phenomenon is clearly visible. 
The voltage bias is not sufficiently good to eliminate the 
leakage, i.e. the P(E)-curvc 8 for tunneling events is not 
sufficiently peaked at the origin. 

These considerations suggest that it is interesting to 
study the difference in the currents AI with the gate 
shifted by 180 degrees. Figure 0Ja) shows the measured 
behavior of AI at 2.5 MHz versus the high level of gate 
voltage with the low level set to zero. The current may 
be seen to increase in clear steps. The expected height 
of a step is twice the pumped current, i.e. 4e/ which in 
this case is some 1.6 pA. Since we sweep the high level of 
the gate signal and not just the amplitude with constant 
offset, the steps should occur at 2e intervals in the gate 
charge. However, due to random parity changes (quasi- 
particle "poisoning" ) at time scales that are much shorter 
than our measurement time scale (0.1 s) but longer than 
the pumping cycle (10 -6 s) we observe the time average 
of two 2e-periodic staircases that are shifted by e in the 
gate charge. For instance in Ref. the tunneling time 
for quasiparticles was estimated to be 10 /is in a similar 
structure while in Ref. ^3 it was some 10~ 2 s for a cou- 
pled system of two superconducting transistors with one 
grounded. We were unable to measure the corresponding 
time in our setup, but based on this supporting evidence 
we argue that the transport of current is due to Cooper 
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FIG. 4: (a) Difference AI in current of forward and backward 
pumping at 2.5 MHz against the high level of the gate signal 
Vg with the low level at zero. The dashed lines are drawn 
at 2e/ intervals, (b) Large gate amplitude behavior of AI at 
a few frequencies. The dashed lines show the expected gate 
dependece, i.e. their slope is 2e/. The curves are offset for 
clarity, (c) Fitted slopes to the data of the previous plots up 
to VgCg/e = 10 are shown by circles. The solid line indicates 
the expected behavior. The voltage bias point was around 10 
\xV in all the above plots. 

pairs since the order in which the -Ej's are manipulated 
changes the direction of current. The quasiparticles effec- 
tively shift the gate charge by e but rarely enough such 
that the pumping is undisturbed on the level of preci- 
sion of the present measurement. If this interpretation 
is made then one sees that the obtained results are in 
very good agreement with theory. Figure E[b) illustrates 
the measured large amplitude behavior of the pumped 
current at frequencies between 1 MHz and 4 MHz. We 
see that the current lacks behind the prediction with in- 
creasing frequency and amplitude. At 1 MHz no clear 



bending of the curve is seen up to gate amplitude of 40e, 
while at 4 MHz the performance starts to degrade after 
lOe. One can observe by looking at Fig. Efb) that the 
" bending" is more pronounced at larger bias voltage val- 
ues (voltage is on the order of 10 /iV in Fig. while no 
visible bending happens up to amplitudes of 68e when 
V ~ 0. Small amplitude behavior in Fig. 01 however, 
is linear aside from the steps with a slope of 2ef. Fig- 
ure 0Jc) shows the slopes obtained from linear fits to the 
data of Fig.0{a) and the ten first steps of Fig.^b). One 
sees that the agreement is again good. 

The above results prove that the flux and voltage 
driven pumping of Cooper pairs is experimentally pos- 
sible in a single-island device. However, in order to serve 
as a practical device the leakage current needs to be taken 
care of as well as the quasiparticle poisoning. The quasi- 
particles may possibly be handled by either quasiparti- 
cle "traps" or by BCS gap profile engineering^. As to 
the reduction of the leakage, several options exist. One 
option is the engineering of the electromagnetic environ- 
ment such that the voltage biasing is good also at fre- 
quencies on the order of the charging energy. This would 
result in DC IV-characteristics heavily peaked at zero 
voltage with negligible leakage current. Another way to 
cut down the leakage is to fabricate a longer chain of 
junctions. A multiloop SQUID would possibly improve 
the suppression of Ej without increasing the number of 
controls. Improved RF-engineering would also be of ben- 
efit in arranging the flux pulses. To conclude, the results 
are encouraging in spite of several nonidealities observed 
and the pumping of Cooper pairs with flux control looks 
much more attractive than with a mere multiple gate 
voltage control. 
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